Percolation on interacting, antagonistic networks 
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Recently, new results on percolation of interdependent networks have shown that the percolation 
transition can be first order. In this paper we show that, when considering antagonistic interactions 
between interacting networks, the percolation process might present a bistability of the equilibrium 
solution. To this end, we introduce antagonistic interactions for which the functionality, or activity, 
of a node in a network is incompatible with the functionality, of the linked nodes in the other 
interacting networks. In particular, we study the percolation transition in two interacting networks 
with purely antagonistic interaction and different topology. For two antagonistic Poisson networks 
of different average degree we found a large region in the phase diagram in which there is a bistability 
of the steady state solutions of the percolation process, i.e. we can find that either one of the two 
networks might percolate. For two antagonistic scale-free networks we found that there is a region 
in the phase diagram in which, despite the antagonistic interactions, both networks are percolating. 
Finally we characterize the rich phase diagram of the percolation problems on two antagonistic 
networks, the first one of the two being a Poisson network and the second one being a scale-free 
network. 
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Over the last ten years percolation proeesses, and more 
in general, dynamical processes in complex netviforks 
P, have gathered great attention. In this context 
it has been shown that complex topologies strongly af- 
fect the dynamics occurring in networks. However, many 
complex systems involve interdependencies between dif- 
ferent networks, and accounting for these interactions is 
crucial in economic markets, interrelated technological 
and infrastructure systems, social networks, diseases dy- 
namics, and human physiology. Recently, important new 
advances have been made in the characterization of per- 
colation [3l-[l0| and other dynamical processes on 
interacting and interdependent networks. In these sys- 
tems, one network function depends on the operational 
level of other networks. A failure in one network could 
trigger failure avalanches in the other interdependent net- 
work, resulting in the increased fragility of the interde- 
pendent system. In fact, it has been shown 0-0] that two 
interdependent networks are more fragile than a single 
network and that the percolation transitions in interde- 
pendent networks can be first order. These results have 
subsequently been extended to multiple interdependent 
networks d, Q and to networks in which only a fraction 
of the nodes are interdependeirt [loj . 

Here, we want to investigate the role of antagonistic in- 
teractions in the percolation transition between interact- 
ing networks. As it is happening in spin systems, where 
antifcrromagnetic interactions can result in the frustra- 
tion of the system, also in interacting network, the pres- 
ence of antagonistic interactions between the nodes intro- 
duce further complexity in the percolation problem. As a 
first step in investigating this complexity in this paper we 
will consider two interacting networks with purely antag- 
onistic interactions. We will show that for two Poisson 
networks with exclusively antagonistic interactions the 
steady state of the percolation dynamics corresponds to 
the percolation of one of the two networks. Neverthe- 



less, the solution of the model is surprising because there 
is a wide region of the phase space in which there is a 
bistability of the percolation process: depending on the 
detail of the percolation dynamics either one of the two 
networks might end up to be percolating. Therefore, in 
this new percolation problem, not only the percolation 
transitions might be first order, but we found that there 
is a real hysteresis in the system as we modify the aver- 
age degrees of the two networks. Furthermore, we extend 
the analysis to networks with other topologies, studying 
the percolation transition in two antagonistic scale-free 
networks, and in two networks one of which is a Pois- 
son network, and the other one is a scale-free network. 
We characterize the rich phase digram of the percolation 
transition in these networks that display both first and 
second order transitions, and, in top of that, might show 
a bistability of the percolation solutions. Interestingly, in 
the percolation phase diagram of these interacting net- 
works there is a region in which both networks percolate 
on the same time, demonstrating a strong interplay be- 
tween the the percolation dynamics and the topology of 
the network. Finally, these results shed new light on the 
complexity that the percolation process acquires, when 
considering percolation on interdependent, antagonistic 
networks. 

Percolation on antagonistic networks. In this paper 
we introduce antagonistic interactions in percolation dy- 
namics on interdependent networks. The difference with 
respect to the case of interdependent networks is that if 
a node i is active on one network it cannot be active in 
the other one. We consider two networks of N nodes. 
We call the networks, network A and network B with 
degree distribution respectively. Each node 

i is represented in both networks. In particular, each 
node has a set of neighbor nodes j in network A, i.e. 
j € N^{i) and a set of neighbor nodes j in network B, 
i.e. j G N^{i). 
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A node i belongs to the percolation cluster of net- 
work A, if it has at least one neighbor j £ N^{i) in the 
percolating cluster of network A, and has no neighbors 
j G N-^{i) in network B that belong to the percolating 
cluster of network B. Similarly, A node i belongs to the 
percolation cluster of network B, if it has at least one 
neighbor j G (i) in the percolating cluster of network 
B, and has no neighbors j G N^{i) in network A that 
belong to the percolating cluster of network A. 

If we define 5*^ as the probability to find a node in the 
percolation cluster of network A, and S's as the proba- 
bility to find a node in the percolation cluster of network 
B, we have 



Sa = [l-G^{l-S'^)]G^il-S's) 



(1) 



where 5*^ is the probability that, following a link in net- 
work A, we find a node in the percolation cluster of net- 
work A, and 5*^ is the probability that, following a link 
in network B, we reach a node in the percolation cluster 
of network B. Moreover, in Eq. ([T]) we have used Gq^^ {z) 

and Gi (z) to indicate the generating functions of net- 
work A and B defined according to the definition 
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where we use the degree distributions p"^ (k) , (k) , re- 
spectively, for network A and network B. According to 
percolation theory, applied in this case to two antago- 
nistic networks, the variables S*^, S'g satisfy, on a locally 
tree like network, the following recursive equations 



(1 - (1 - S'^))G^{1 - S's) = fAiS'A, S's), 



S'b = {l-Gf{l-S'B))G^{l-S'A)^ fB{S'A,S's).{3) 

The solutions to the recursive Eqs. ^ can be classified 
into three categories: 

• (i) The trivial solution in which neither of the net- 
work is percolating 5^ ~ S'g ~ 0. 

• (ii) The solutions in which just one network is per- 
colating. In this case we have either 5^ > 0, S*^ = 
or 5^ > 0, 5"^ = 0. From Eqs. ^ we find that the 
solution S'j^ > 0, S'g = emerges at a critical line 
of second order phase transition, characterized by 
the condition 
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(4) 



Similarly the solution S'^ > 0,5^ = emerges 
at a second order phase transition when we have 

^^"''(fc>fl ''^ ~ ^' "^^^^ condition is equivalent to the 
critical condition for percolation in single networks, 
as it should, because one of the two networks is not 
percolating. 



(Hi) The solutions for which both networks are per- 
colating. In this case we have > 0, 5*^ > 0. 
This solution can either emerge (a) at a critical line 
indicating a continuous phase transition or (b) at 
a critical line indicating discontinuous phase tran- 
sition. For situation (a) the critical line can be 
determined by imposing, for example, S*^ — > in 
Eqs. ([3]), which yields 



5' = 



l-Gf(l-5^), 
{k{k-l))A 



G^{l-S's). 



(5) 



A similar system of equation can be found by using 
Eqs. ([3]) and imposing S'^ — > 0. For situation (b) 
the critical line can be determined imposing that 
the curves S'^ = fA{S'A, S'g) and S's = fniS'^, S'^), 
are tangent to each other at the point where they 
intercept. This condition can be written as 



dfA 
dS'A 



fdfE 



\dS' 



dfA dfB 
dS'g dS'^ 



0, 



(6) 



where 5^, Sg must satisfy the Eqs. ([3]). 



The stability of the solutions. Not every solution of the 
recursive Eqs. ([3]) is stable. Therefore, we check the sta- 
bility of the fixed points solutions of Eqs. by lineariz- 
ing the equations around each solution. The Jacobian 
matrix J of the system of Eqs. ([3]) is given by 



J = 



dfA dfA 

ds'. as' 

dfB dfB 

dS'. dS' 



(7) 



The eigenvalues Ai^2 of the Jacobian can be found by 
solving the characteristic equation | J — A/| ~ 0, which 
reads for our specific problem. 



dfA 
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dfA dfB 
dS's dS'^ 



= 0. 



(8) 



A solution is stable if and only if Ai_2 < 1- Assuming that 
the eigenvalues of the Jacobian corresponding to each 
solution of the Eqs. ([3]) change continuously when we 
smoothly change the parameters determining the topol- 
ogy of the networks, the change of stability of each solu- 
tion will occur when max(Ai , A2) = 1. In the following we 
will discuss the stability of the solutions of type (i)-(iii). 

• (i) Stability of the trivial solution S'j^ = S'^ = 0. 
The solution is stable as long as the following two 
conditions are satisfied: 



A1.2 — 



{k{k-l))A/l 
{k)A/B 



< 1. 



(9) 



Therefore the stability of this solution change on 



the critical lines 
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FIG. 1: Phase diagram of the percolation process on two an- 
tagonistic Poisson networks of average degree {k)A ~ za and 
(fc) B ~ zb respectively. 



• (a) Stability of the solutions in which only one net- 
work is percolating. For the case of 5^ 
the stability condition reads 



OS'b>0 



Ai 
A2 



Gf(z) 



dz 



< 1 



--1-S' 



(10) 



We note here that if A2 > Ai we expect to observe a 
change in the stability of the solution on the critical 
line given by Eqs. ([5]). A similar condition holds 
for the stability of the solution S*^ > 0,Sg ~ 0. 

• (Hi) Stability of the solution in which both networks 
are percolating 5^ > 0, 5*^ > For characterizing 
the stability of the solutions of type (iii) we have to 
solve Eq. (jH) and impose that the eigenvalues Ai^2 
are less then 1, i.e. Ai^2 < 1- We observe here that 
for A = 1 Eq. ^ reduces to Eq. Therefore we 
expect to have a stability change of these solutions 
on the critical line given by Eq. ([6|). 

Two Poisson networks. In order to consider a specific 
example of antagonistic networks we consider two antag- 
onistic Poisson networks with average degree (fc)^ = 
and (fc)_B = Zb respectively. In the following we charac- 
terize the phase diagram of the percolation dynamics on 
two antagonistic Poisson networks shown in Fig. [T] 

In region I, with za < zb < 1, the only stable solu- 
tion is the trivial solution S'^ — S'g = 0. 

In regions II-A, with za > 1 and 1 < zb < ln(2:^) /(I — 
1/za), the only stable solution is S'a > 0,S'g = 0. Net- 
work A has a larger average connectivity than network B 
and is percolating, while network B is non percolating. 

In regions II-B, with > 1 and 1 < za < ln(zB)/(l — 
1/zb), the only stable solution is S'g > 0,S'^ = 0. Net- 
work B has a larger average connectivity than network A 
and is percolating, while network A is non percolating. 
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FIG. 2: (Color online) Panels (a) and (b) show the hysteresis 
loop for the percolation problem on two antagonistic Poisson 
networks with zb = 1.5. Panels (c) and (d) show the hys- 
teresis loop for the percolation problem on two antagonistic 
networks of different topology: a Poisson network of average 
degree za = 1.8 and a scale-free networks with power-law 
exponent 7b, minimal degree m = 1 and maximal degree 
K = 100. The hysteresis loop is performed using the method 
explained in the main text. The value of the parameter e used 
in this figure is e = 10"'^. 



In region III, with za > ln(zB)/(l — 1/zb) and zb > 
ln(zA)/(l — ^/za) , there are two stable solutions S'^ > 
0,S'b = and S'g > 0,S'^ = 0. In this region, both 
network A and network B might be percolating. The 
system is bistable and Eqs. (jS]) has two possible solutions. 

One should note that the solution S'a > 0, S'^ > in 
which both networks are percolating is always unstable 
in this case. This implies that at any given time only one 
of the two networks is percolating. 

In order to demonstrate the bistability of the percola- 
tion solution in region III we solved recursively the Eqs. 
([3]) for = 1.5 and variable values of za (see Figured]). 
We start from values of za = 4, and we solve recursively 
the Eqs. ©. We find the solutions S'a = S'a{za = 4) > 0, 
S'g — S'g{zA = 4) = O.Then we lower slightly za and we 
solve again the Eqs. ([3]) recursively, starting from the ini- 
tial condition S'X = S'^izA = 4) -he, S'£ = S'g{zA = 4) -He, 
and plot the result. (The small perturbation e > is 
necessary in order not to end up with the trivial solution 
S'ji^ = 0,5^ = 0.) Using this procedure we show that if 
we first lower the value of za and then again we raise it, 
spanning the region III of the phase diagram as shown 
in Figure [21 the solution present an hysteresis loop. This 
means that in the region III either network A or net- 
work B might end up to be percolating depending on the 
details of the percolation dynamics. 

Two scale-free networks. Here, we characterize the 
phase digram of two antagonistic scale- free networks with 
power-law exponents 7a, 7b, as shown in Figure [31 The 
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FIG. 3: The phase diagram of the percolation process in two 
antagonistic scale-free networks with power-law exponents 
7a, 7s for finite networks. The minimal degree of the two 
networks is m = 1 and the maximal degree K. The panel on 
the left show the effective phase diagram with K = 100 and 
the panel on the right show the effective phase diagram for 
K = 10,000. 



two networks have minimal connectivity m = 1 and vary- 
ing value of the maximal degree K. 

The critical lines of the phase diagram depend on the 
value of the maximal degree K of the networks. There- 
fore, the plot in Figure [3] has to be considered as the 
effective phase diagrams of the percolation problem on 
antagonistic networks with a finite cutoff K. The phase 
diagram is rich, showing a region (Region III) in the fig- 
ure where both networks are percolating demonstrating 
an interesting interplay between the percolation dynam- 
ics and the topology of the network. 

A description of the phase diagram follows. 



In region I, with \k)l ^ ^ ^^"^ 
trivial solution S*^ = 0, 5^ = is stable. 

In regions II-A the solution in which network A perco- 
lates S'j^> 0,S'g = is stable. 

In regions II-B the solution in which network B perco- 
lates S'^ > 0, 5^ = is stable. 

In region III the solution in which both network A and 
B percolate S'^ > 0, S'^ > is stable. 

A Poisson network and a scale-free network. Finally 
we consider the case of a Poisson network (network A) 
with average connectivity {k)A = za, and a network B 
with scale-free degree distribution and power-law expo- 
nent of the degree distribution 73. The scale- free net- 
work has minimal connectivity m = 1 and maximal de- 
gree given by K . In Figure |4] we show the phase diagram 
of the model in the plane {-^BtZa)- The critical lines 
of the phase diagram are dependent on the value of the 
cutoff K of the scale-degree distribution. For these rea- 
sons we have to consider the phase diagrams in Figure 
m as effective phase diagrams of the percolation prob- 
lem on networks with maximal degree K. The phase 
diagram includes two regions, (region III and region V) 
with bistability of the solutions and two regions (region 
IV and region V) in which the solution in which both 
networks are percolating is stable. In the following we 
describe the percolation stable solution in the different 
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FIG. 4: (Color online) Phase diagram of the percolation pro- 
cess on a Poisson network with average degree (fc)A = za 
interacting with a scale-free network of power-law exponent 
7_B, minimal degree m = 1 and maximal degree K. The panel 
on the left show the effective phase diagram for K = 100 and 
the panel on the right show the effective phase diagram for 
K = 10,000. 



regions of the phase diagram in Figure 21 

In region I, with ZA<1, ''''%)]^^'' < 1, the only stable 
solution is the trivial solution S'y^ = S'^ = . 

In region II-A the solution for which network A is per- 
colating 5^ > 0, 5^ = is stable. 

In regions II-B the solution for which network B is 
percolating S*^ > 0, S*^ = is stable. 

In region III there are two stable solutions, in which 
cither network A percolates S"^ > 0, Sg = or network 
B percolates S'^ > 0, 5^ = 0. Therefore in this region 
we have bistability of the solutions. 

In region IV there is only one stable solution in which 
both networks percolate, 5^ > 0, S'^ > 0. 

In region V there are two stable solutions. In the first 
one, network B percolates S'g > 0, 5^ = 0, and in the 
second one, both networks percolate S'j^ > 0, S'g > 0. In 
this region there is a bistability of the solutions of the 
percolation problem. 

In order to demonstrate the bistability of the perco- 
lation problem we solved recursively the Eqs. ([3]) for 
zb = 1.8 (see Figure[2]). We start from values of 75 = 3, 
and we solve the Eqs. We find the solutions S'^ = 
S'a{ib = 3) > 0, S"^ = S'^(7B = 3) = 0. Then we lower 
slightly "fB and we solve again the Eqs. ^ recursively, 
starting from the initial condition S'^ = S'^{'Jb = 3) + e, 
— '^'si^B = 3) + e and we plot the result. (The small 
perturbation e > is necessary in order not to end up 
with the trivial solution S'^ = 0,5^ = 0.) Using this 
procedure we show that if we first lower the value of 73 
and then again we raise it as shown in Figure [H the so- 
lution present first a second order phase transition to a 
phase in which both networks are percolating and then 
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an hysteresis loop in correspondence of region V. This 
demonstrates the bistabihty of the solutions in region V 
and the existence of a phase in which both network per- 
colate in region IV and region V. 

Conclusions. In conclusions, we have investigated how 
much antagonistic interactions modify the phase diagram 
of the percolation transition. The percolation process on 
two antagonistic networks shows important new physics 
of the percolation problem. In fact, the percolation pro- 
cess in this case shows a bistability of the solutions. This 
implies that depending on the details of the percolation 
dynamics, the steady state of the system might change. 
In particular, we have demonstrated the bistability of the 
percolation solution for the percolation problem on two 
antagonistic Poisson networks, or two antagonistic net- 
works with different topology: a Poisson network and a 



scale-free network. Moreover, in the percolation transi- 
tion between two scale-free antagonistic networks and in 
the percolation transition between two antagonistic net- 
works with a Poisson network and a scale-free networks, 
we found a region in the phase diagram in which both 
networks are percolating, despite the presence of antago- 
nistic interactions. We believe that this paper opens new 
perspectives in the percolation problem on interdepen- 
dent networks, which might include both interdependen- 
cies and antagonistic interactions eventually combined in 
a boolean rule. In an increasingly interconnected world, 
understanding how much these different types of inter- 
actions affect percolation transition is becoming key to 
answering fundamental questions about the robustness of 
interdependent networks. 
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